Correlation between stick-slip frictional sliding and charge transfer 
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A decade ago, Budakian and Putterman (Phys. Rev. Lett., 85, 1000 (2000)) ascribed friction to 
the formation of bonds arising from contact charging when a gold tip of a surface force apparatus 
was dragged on polymethylmethacrylate surface. We propose a stick-slip model that captures the 
observed correlation between stick-slip events and charge transfer, and the lack of dependence of the 
scale factor connecting the force jumps and charge transfer on normal load. Here, stick-slip dynamics 
arises as a competition between the visco-elastic and plastic deformation time scales and that due 
to the pull speed with contact charging playing a minor role. Our model provides an alternate basis 
for explaining most experimental results without ascribing friction to contact charging. 

PACS numbers: 81.40.Pq, 46.55.+d, 05.10.-a 
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I. INTRODUCTION 

Despite its long history, several aspects of friction re- 
main ill understood even today. This can be partly at- 
tributed to the fact that mechanisms contributing to fric- 
tion are scale dependent and those operating at lower 
scales contribute to friction at higher scales. The situa- 
tion is further complicated as several other factors such 
as the possible presence of interfacial layer between the 
contacting surfaces, plastic deformation of the contacting 
asperities, contact electrification etc., are also known to 
contribute. Recent advances in friction force microscope 
and surface force apparatus (SFA) have provided useful 
tools to understand friction at these scales [iHj]. This 
coupled with large scale molecular dynamics (MD) sim- 
ulations [g, |6| have provided a better understanding of 
the mechanisms that contribute to friction at nanometer 
contacts and their relevance to macroscopic friction. Ex- 
perimental studies show that Amontons' law (the linear 
relationship between frictional force and load) holds in a 
number of situations (from nanometer to macroscopic di- 
mensions) when adhesion is small or the surface is rough 
(either dynamically induced or otherwise [l|-|l]). How- 
ever, smooth adhesive contacts follow Johnson-Kendal- 
Roberts (JKR) theory [7], Q. Recent MD simulations 
support experimental findings to a large extent [5|, |9| . 

In contrast, the dynamical aspects of intermittent fric- 
tional sliding are less studied 0, S 03 ■ Efforts to under- 
stand stick-slip dynamics have established that plastic 
deformation of the contacting interfacial material is re- 
sponsible for slip |Q|— fl6Jj . A decade ago, Budakian and 
Putterman [17 ] , proposed a new point of view that fric- 
tion is due to the formation of bonds arising from contact 
charging. This claim however is at variance with ear- 
lier studies that show n egligibl e contribution from con- 
tact charging to friction |18l - [21| for charge density ~ 10 8 
charges/mm 2 reported in Ref. 



17] 



The authors demonstrate a correlation between stick- 
slip events and charge transfer when the gold tip (ra- 
dius 0.5 mm) of the cantilever of a SFA was dragged 
with a velocity of a few fim/s on polymethylmethacrylate 
(PMMA) surface [17[. Typical measured charge density 



~ 10 8 charges/mm 2 . The magnitude of the slip events is 
proportional to the ensuing charge transfer to the PMMA 
surface. Interestingly, the total force and the total charge 
deposited over a scan length collapses onto a single curve 
when an appropriate choice of a scale factor a is made. 

o 

The value of a ~ 0.4 eV/A is close to the energy win- 
dow for transfer of charge between the surface states of 
PMMA and metallic Fermi level (H). Further, the au- 
thors find a to be nearly constant in the range of normal 
loads F n from 68 to 106 mN. Based on these obser- 
vations, they ascribe frictional force to the formation of 
bonds due to contact charging at the interface. Snapping 
of bonds leads to slip. To the best of our knowledge, there 
is no model that explains these intriguing results. Our 
purpose is to construct a model that captures the major 
results [17]. Our model provides an alternate explanation 
for the observed correlation and the lack of dependence 
of the scale factor a on normal load. The threshold for 
slip is determined by plastic deformation of the interface 
material [id 1 121 ] while contact charge plays a minor role. 



II. BACKGROUND 

A. Motivation and approach 

A notable feature of stick-slip systems is that the stick 
phase lasts much longer than the slip phase. Then, it 
is intuitively clear that charge builds up during the stick 
phase and charge transfer to the PMMA substrate occurs 
during the slip phase. In view of this, our idea is to first 
construct a stick-slip model based on contact mechan- 
ics of single asperity contact and other relevant physical 
processes, and couple it to charging and charge transfer 
equations. Then, as the contribution from contact charg- 
ing to friction is small [l8l - t23 | , we expect the correlation 
between slip events and charge transfer should follow. 
Indeed, if one is interested in bringing-out the correla- 
tion between stick-slip events and charge transfer only, a 
'toy' model for stick-slip that includes velocity weakening 
law as an input coupled to charging and charge transfer 
equations is adequate. Such a model has been shown 



to reproduce the required correlation [23|. However, as 
our interest is to recover most experimental features, we 
attempt to include all relevant features of friction and 
contact mechanics. 

Sliding friction is well recognized to be a complex phe- 
nomenon with multiple possible scenarios depending on 
the precise experimental conditions [lj, |9| . The nature of 
the results depend on variety of factors such as the size 
of the tip radius (typically 10 — lOOnrn radius for AFM 
and mm — cm radius for SFA), the load levels, the nature 
of contacting materials, the nature of the initial and fi- 
nal microstructures of the tip and the substrate (smooth 
or rough) , the magnitude of adhesion or the presence of 
intervening layer etc[l|-|4J,|9|]. For example, very different 
conditions such as low adhesion and rough or damaged 
surfaces can lead to similar friction- load relationship 1- 
4]. For our case, such detailed information on the precise 
experimental conditions on the microstructure of the gold 
tip or the nature of surface of the PMMA substrate be- 
fore and after the scan, is not provided in Ref. [l7|- In 
view of this and in view of the fact that any model build- 
ing effort to understand the stick-slip dynamics requires 
quantitative information on the contact area, we use con- 
tact mechanics as the basis of our model. Moreover, the 
JKR theory has been validated in a number of SFA like 
situations. 



B. Single asperity contact, visco-elasticity, 
plasticity and contact charging 

By necessity our approach to modeling will be dynam- 
ical. The idea is to construct a stick-slip model that 
builds on known results for elastic deformation of the 
contacting surfaces by adding contributions from visco- 
elastic and plastic deformations of the contacting sur- 
faces. Contact mechanics offers a mathematical basis for 
quantitative description of contact radius and penetra- 
tion depth. The theory deals with smooth nonadhesive 
and adhesive contacts is designed as an equilibrium the- 
ory. Yet, experiments show that it is applicable to a num- 
ber of sliding conditions as wellfll la [9. 124] ■ Non-adhesive 
contact was addressed by Hertz later extended to adhe- 
sive contacts by Johnson-Kendal-Roberts and later by 
Maugis 0, S 123- In these theories, the area of the con- 
tact is a sub-linear function of the normal load F n |7I.|25|. 
As the Hertz contact area differs less than 5% from the 
JKR contact area (for normal loads F n and tip radius R 
used in [17|), we use the Hertz's contact area A n given 
by An = 7to 2 = Tr[3RF n /4E*} 2 / 3 = nRz. Here, a is the 
contact radius, z is the penetration depth, and E* is the 
effective elastic constant (of the contacting materials). 
The value of a is ~ 23/im for load F n — 0.1N. Thus, the 
estimated value of stress is ~ 60MPa that is higher than 
the compressive yield stress of PMMA [2q |. but smaller 
than the yield stress of gold ~ 80MPa. (See Table 1.) 
Thus, the plastic deformation is largely confined to the 
weaker PMMA material. This view is supported by MD 



simulations as well [9|, iRfl - flq ] . 

Quite early, it was shown that accelerated creep of the 
asperities leads to slip [lj| • Later studies show that stick- 
slip instabilities arise from two competing mechanisms, 
namely, 'geometric aging' of the contacting asperities and 
'rejuvenation by motion due to plastic deformation' [111 ]. 
Models for stick-slip instability have been proposed using 
these ideas [H| LUJ- Indeed, aging kinetics of contacting 
asperities has been used to derive a 'N'-shaped veloc- 
ity dependent friction coefficien t Il2l . a generic feature of 
most stick-slip systems [lCj, [llj, |27| . Dislocation assisted 
model for frictional sliding has also been suggested [281 ]. 
Further, MD simulations also show plastic deformation 
of the contacting interfacial material and possible trans- 
fer of material between the contacting surfaces [9l 1 1314 la ] . 
Thus, we consider plastic deformation of the interfacial 
layer is responsible for slip. 

Even though the subject of contact electrification is 
rather old, the mechanisms underlying the charge trans- 
fer between a metal and polymer are still debated [2lll23 |. 
However, the fact that contact electrification contributes 
to adhesion is well recognized. A number of early stud- 
ies suggest small contribution to adhesion [l8l - [22l ] for the 
charge levels reported in Ref. [13] • An estimate of the 
adhesive force from contact charging can be obtained by 
considering a charged double layer formed at the inter- 
face. Noting that contact charging can occur only at the 



area of contact given by A r , 



ttRz, the attractive 



force is given by irRza 2 /2eoK, where a is the charge den- 
sity and eo is the vacuum permittivity and k the dielectric 
constant. Using the contact radius (a ~ 23 fim for nor- 
mal load O.liV), the force of attraction is ~ 10 _9 iV that is 
several orders of magnitude smaller than mN force drops 
observed in experiments [17] . 

Here, it is pertinent to point out that the bond forma- 
tion attributed to contact charging is very different from 
conventional bonds formed when two atomically 'clean' 
flat surfaces are brought together or when two surfaces 
are pressed together. The bonds so formed would have 
all the relevant electronic contributions (such as kinetic 
energy of the electron charge density, the long range elec- 
trostatic interaction, exchange energy, correlation energy 
etc). In addition, there could be charge transfer also when 
the Fermi energies of the two metals are different or when 
a metal is in contact with a suitable polymer, as in this 
case. Indeed, the process of 'bond formation' is well mim- 
icked by MD simulations that use appropriate 'potentials' 
by switching-on the interaction between atoms of the two 
surfaces when they are brought into contact. 

To summarize, we list the basic ingredients of our 
model. First, as stated earlier, we assume that the 
mechanics of single smooth asperity in contact with a 
smooth surface is valid for the current situation 0, H|. 
Second, we note that PMMA is a much softer material 
compared to gold and is a visco-elastic material so that 
we include visco-elastic contribution for the deformation. 
Third, as shown above, for the load levels and tip radius 
in the experiment, stress level can exceed the yield stress 



of PMMA, but remains less than that of gold. Thus, 
we assume that the dissipation is confined to the weaker 
PMMA material. (Note however, due to the mean field 
nature of our model, our model does not have any scope 
for dealing such details as where the dissipation occurs.) 
Finally, we include the frictional resistance arising from 
contact electrification. Thus, in view of the fact that 
stick-slip dynamics arises from the interplay of all inter- 
nal relaxational time scales with the applied time scales, 
we expect that the inclusion of time dependent contribu- 
tions arising from visco-elasticity and plastic deformation 
of the contacting interface material will lead to stick-slip 
dynamics. 



III. MODEL EQUATIONS 

Our stick-slip model (similar to that of Ref. 29] ) de- 
scribes the center of the contact area x and penetra- 
tion depth z. Our equations are of the general form 
my = F a — F r , where y represents x or z, m is the mass of 
the gold tip, F a is the applied force (shear force F or nor- 
mal force F n ), and F r is the total material response. In 
equilibrium, F a = F r . However, in dynamic conditions, 
there is always an imbalance between the applied force 
and material response due to time dependent responses 
from visco-elastic and plastic deformation. 

Consider motion in the x direction. The reaction force 
F r is the sum of the frictional force Ff, the adhesive force 
from contact charging, and time dependent contributions 
from visco-elastic and plastic deformation. Frictional 
force exerted by the tip on the substrate is Ff = tqA x , 
where tq is interfacial shear strength and A x is the pro- 
jected area in the a;— direction. The latter is obtained by 
assuming a parabolic tip (defined by x 2 — 2Rz) that gives 
A x ~ ^nV^Rz 3 / 2 = A X fiZ 3 / 2 . In our model, we include 
both visco-elastic and plastic deformations of interfacial 
material. Visco-elastic effects are generally included by 
replacing the elastic stress r by time dependent visco- 
elastic stress r + ye, where r] is the viscosity and e is the 
strain rate [29|, [3(j. Then, the visco-elastic creep con- 
tribution is mA x x/D — rinA Xt oz 3 ' 2 x/D, where rjn is the 
shear viscosity in the a;— direction with D a length scale 
to be fixed. 

The contacting asperities undergo plastic deformation 
beyond the yield stress r y . The general phenomenologi- 
cal expression for plastic strain rate of crystalline mate- 
rials is given by e p = eoexp[— J^ Va ] , where AF is the 



kT 

change in the free energy, V a is the activation volume, 
r is the stress, and eo is an appropriate prefactor. This 
expression has been adopted for the case of polymers as 
well 26]. However, this contains several parameters that 
are unknown for PMMA. An alternate phenomenologi- 
cal expression for plastic strain rate of crystalline ma- 
terials that has been used is given by e p = eo(r/r a ) n , 
where r is stress, r y is the yield stress of the mate- 
rial, n an exponent, and eo is the strain rate at some 
specified value of the stress [23, |3IH33| . This expres- 



sion has been successfully used even in the case unsta- 
ble intermittent plastic flow [23, [32|, [33[- To minimize 
the number of parameters, we have adopted this expres- 
sion for plastic deformation of the interface material. In 
the physical context, we identify the threshold stress To 
with the yield stress of the softer PMMA material. Not- 
ing that the deformation is continuous beyond the linear 
visco-elastic flow into the nonlinear plastic flow regime, 
the total contribution from these two flows can be writ- 
1 - (- — ^ttw) 1 EH. It is clear 
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A X fiTQZ 3 / 2 , the flow is resistive, 
changes over to strain rate softening flow when F ex- 
ceeds the threshold force A x t z 3 / 2 as the flow rate in- 
creases abruptly. Finally, we include the frictional resis- 
tance from contact charging given by /iA n a 2 /2eqk with 
fi referring to the frictional coefficient. We note here that 
the magnitude of the contact charge contribution (which 
is ~ nN for the charge level in [17j ) is several orders of 
magnitude less than the magnitude of the force drops ( 
of mN ) in experiments. In contrast, the frictional re- 
sistance that is ~ lmN as can be verified by using the 
value of r — lOMPa and z = a 2 / R in A xfi T z 3 / 2 (Ta- 
ble HJ. However, the frictional resistance from contact 
charging can become important when the contact charge 
density is high. Then, the inclusion of this contribution 
increases the threshold force for slip from A x ^tqz 3 / 2 to 
A Xj qTqz 3 ' 2 + unRza 2 /2k€q. Then, the equations for x 
and F are 
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A x ,oToz 3 / 2 + nirRza 2 /2k,6q 



F 



K\\(V a -x), 



(2) 



where A(z) = A xfi r z 3 / 2 . In Eq. Q, F = -K\\ (x - V a t) 
is the applied force with t, K\\ and V a referring respec- 
tively to time, the effective lateral spring constant and 
applied velocity. Eq. ^ is the differential form of F. 
We further make two choices. The first choice is to use 



A(z) = Aq, a constant 



Most of the results pre- 



sented here are for Model I. The second choice is to re- 
tain A(z) — A x _ z 3 / 2 . We refer to this as Model II. As 
we shall show the results of Model II are similar to those 
of Model I. 

Following a similar approach, the equation for z is the 
difference between the normal force F n and the upward 
response of the material F r . This is the sum of elas- 
tic contribution and contributions arising from the time 
dependence of visco-elastic and plastic deformation pro- 
cesses. The elastic response is given by ^R}l 2 z 3 l 2 E* . 
Note that in equilibrium, the balance between the normal 
load F n and the elastic response gives the Hertz relation 
a 2 = Rz = {ZRF n /AE*) 2 / 3 . Since, the substrate mate- 
rial is a visco-elastic material, to account for the visco- 
elastic response, we replace E* by E* + %i, where rj± 
refers to the bulk viscosity of the PMMA [30||. Further, 
due to the normal load, there can be plastic deformation 



(of PMMA) when the load and contact area are such 
the stress r = F n /i:Rz is larger than the compressive 
yield stress T y , n . Again using the phenomenological rela- 
tion e p = eofr/Ty) 9 , we can write the plastic strain rate 

contribution to be — ^jz( - — J^ j . q is an appropriate 
exponent for the compressive deformation. Finally, while 
the tip creeps in the z direction, contact area and hence 
the position of tip center x, also creep in the x direction 
due to shear force F. Noting that z = z(x(t),t), we get 
an additional creep contribution crj±x/D. Here J§ = c is 
taken to be a constant. (It must stated that the depen- 
dence of z on x is not known. In general ^ would be a 
function of z which has been assumed to be constant in 
this case.) Combining these contributions, we have 
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The last term is the normal component of the shear force. 
Eqs. (I1I2I3I) support stick-slip dynamics for a range of 
parameter values. 

In our model, stick-slip instability arises due to a feed 
back loop. To see this, we note that when x is small 
(stuck phase), F increases (see Eq. ©). Concomitantly, 
z approaches the equilibrium value at a rate controlled by 
r]±. Thus, the frictional resistance tqAx^z 3 ^ 2 and contact 

2 

charging resistance uttRz-^- increase. Since x < V a , 
F keeps increasing till the threshold for slip is crossed 
triggering an abrupt increase in x due to plastic flow 
(fifth term in Eq. (fT])). This causes the force to drop as 
x overshoots beyond V a during the slip duration (see Eq. 
©). The whole process starts all over again since the 
force has dropped below the threshold for slip. 

Assuming a finite time scale for charging after contact 
is established between the tip and PMMA, the equation 
for the charge density a at the contacting surface is 



a = — I 1 a, 

t„. V a m J D 



(4) 



where a m is the saturation value and t a a time con- 
stant. Once slip occurs, contact charge is transfered to 
the PMMA surface at a local slip rate denoted by the last 
term in Eq. (|4j . Using Eq. (|4|) and the fact that charging 
can only occur at the area of contact determined by the 
load, we get the equation for the total charge a t = irRza, 
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The last term is the charge transferred to the PMMA 
surface during a slip event. Then, the charge transferred 
during each slip event is given by a a = xat/D. 



A. Scaled Equations 

We cast the equations for Model I (i.e., A(z) = Aq) 
in a dimensionless form using a basic length scale D = 
Fmax/Kn and a time scale determined by to 2 = Kn/m. 
Then, the scaled variables are : x = XD, z = ZD. 
F s = F/K^D, fo = toA^D 1 / 2 /^, F n = FJK^D, 
the normal spring constant K± = ^{RD) X / 2 E* , fj± — 
^(RD) 1 / 2 ur) A _/K\\, 7/|| = Aori\\u/K\\D, f y = T y , n irR/K\\, 
v a = Va/tdD, and T a = t a u>. Defining the scaled total 
charge S = ttRDZct/(t with a = {2ttk€ RD 2 K^) 1 / 2 
and S m = irRDa m /ao, the dimensionless equations are: 
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The scaled equation for charge transferred (to the sub- 
strate) is Srf = XY,. As there is only one way coupling 
between (X, Z) variables and E and S^, the instability 
domain of Eqs. (|6I7I8I9[) are nearly the same as the stick- 
slip model (Eqs. ((IPI ISJI) for a m - 10~ 5 C/m 2 reported 
in Ref. [17| . However, when the charge levels are much 
higher, the threshold for slip increases and the instability 
domain is altered from that at low levels of charge. For 
Model II, the scaled equation for X is given by 



X = F s -? Z 3/2 -fi 
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( 10 ) 

where rf« = a x ,oD 1 / 2 ijjr]\\/K\\. Note that except for r)',,, all 
other parameters are the same. Most numerical results 
are for n = q = 1. For the sake of completeness, we have 
presented results for n = 1.75, q = 2 also. 

Table U shows the parameter values (for Model I). Val- 
ues in the second row are from [17J. As no material 
parameters are given in [17] . values of unsealed param- 
eters (fourth row) are taken from the literature wher- 
ever available. In particular, viscosity r\ that depends 
on molecular weight, temperature, shear rates etc, is not 
available. The magnitude of tq ~ /ir s , where r s is the 
shear yield stress and \i ~ 0.3 is the friction coefficient. 
Yield stress depends on strain rate [26| . In addition, (for 
PMMA valid also for polymers in general) compressive 
yiel d stress is always higher than the shear yield stress 
[2a ] . The range of values of these quantities are listed in 
the fourth row. The scaled parameters used for the cal- 
culations are in the last row. Other parameters used are: 
A - a 2 - lCr 10 m, k = 3.3, v a = 1.45 x 1CT 5 , T a = 2 and 



0.122. Choosing F„ 



47 mN gives D = 10" 



and thus the range of F n corresponding to 
is 1.446-2.26. 



106mN 



TABLE I: Parameter values used for the model. Values in the 
second row are from Ref. 11711. See text. 



R 

(mm) 


m 
{kg) 


K \\ 
(N/m) 


(mN) 


V a 

(pm/s) 


(C/m 2 ) 


0.5 


10" 5 


47 


68 - 106 


~ 10 


1.67 x 10" 5 


E* 
GPa 


m 

Pa.s 


Pa.s 


TO 

MPa 


Ty,n 

MPa 


o~o 
(nC/rn 2 ) 


1-3 






0.1 - 10 


1-50 


3.0 


K± 


*?|| 


Vx 


TO 


f y 


^m 


~10 6 


10 2 


8 x 10 8 


1.0 


2500 


0.0167 



Eq. ([6]) (or Eq. (|10l) ) is small compared to the frictional 
term. (See Tableland note that S approaches its satura- 
tion value given by T, m Z asy , where Z asy is the near satu- 
ration value of Z determined by the balance between F n 
and K±Z 3 ' 2 /Kn .). Thus, we can conveniently examine 
the stability of Eqs. ([517191 (or Eqs. (|10I7|9|) ) by drop- 
ping the charge contribution. A standard way is to carry- 
out the stability analysis of these equations in the steady 
state sliding conditions (i.e., X — v a ,X ~ Q,Z = and 
Z = valid in most piezo controlled SEA situations) is 
to find the eigen values of the linearized system of equa- 
tions. However, it turns out that even the calculation 
of fixed points is messy and the expressions are compli- 
cated, forcing us to carry out the analysis numerically 
(not presented here). Instead, a simpler method is find 
F s as a function of the pull velocity [29] in the steady 
state. Setting X = in Eq. © (or [10]), both Model I 
and II give F s = f Q Z 3 / 2 . Setting Z = Z = 0, we get 




Z 3/2 



F„ 



time 




FIG. 1: (color online) (a) Plots of scaled force and cumulative 
charge transferred to the substrate, (b) Collapse of the mean 
displacement X and charge transferred E^ for F n — 2.0 for a 
scale factor a s = 3.37 x 10 . 



IV. RESULTS 

We first examine the stability of the uniform sliding 
state for Model I (Eqs. (|6I7I8I9|) ) or for Model II (Eqs. 
(|10|7|8I9|) ). Here also we use n = q = 1 for simplic- 
ity in calculation. As stated earlier, for charge levels of 
a rs> 10~ 5 C/m 2 , the contribution from the third term in 



[§7 + CVA-Va] 



(11) 



where to simplify the calculations, we have dropped the 
term \ ^F s Z 1 / 2 in Eq. (0 as this term is small com- 
pared to §^Z 3 / 2 . Using this in F s = f Z 3 / 2 , we get 
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r F n 
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Clearly ^f— is negative even from very small values of 
the pull velocity v a . Thus, the steady sliding state is 
unstable. However, there is a lower cut off in v a that can 

be calculated once 

stick-slip oscillations for the systems of equations. Note 

that the static friction threshold is obtained by setting 

X = v a = which gives F s = -^? — -. The analysis also 
shows that the static friction is larger than the sliding 
friction. 

We now consider the numerical solution of Eqs. 
(I6I7I8I9I) . The force-displacement curve is shown in Fig. 
□Ja) for F n = 2.0 (F n = 94 mJV) along with the cumu- 
lative charge transferred to the substrate E^. Clearly, 
the mean force F s gradually increases and saturates, a 
feature seen in experiments (Fig. 3(a) of Ref. [17|)- Fur- 
ther, the average slope of E^ as a function of time (in 
the asymptotic regime) is proportional to that for the 
displacement X. Using a proper scale factor a s these 
two curves can be made to collapse onto a single curve 
shown in Fig. [l](b) for a s = 3.37 x 10 4 . Thus, the model 
recovers the correlation between the stick-slip events and 
charge transfer. 

We now examine the influence of the normal force 
F n . We have calculated displacement X and cumula- 
tive charge transferred E^ for F n = 1.446 — 2.26 (F„ = 
68 — 106 mN). In each case, X and E^ curves collapse 
onto a single curve for a proper choice of the scale factor 



a s . In particular, for F n = 1.446 and 2.26 the values of 
a s respectively are 4.18 x 10 4 and 3.09 x 10 4 . Then, the 
ratio of the maximum deviation (1.09 x 10 4 ) to the mean 
(3.635 x 10 4 ), is ~ 30%. This compares well with 25% 
scatter in a in experiments for F n = 68 — 106 mN (Fig. 
3(b) inset of [l7[). This shows that the change in a s 
(equivalently a) is small for the limited range of normal 
loads considered (F n = 1.446-2.26 or F n = 68-106miV). 

In our model, as A n ex F n , we expect the scale factor 
a to depend on the load. However, the area changes by a 
factor 1.34 when F n is changed from 68 to 106 mN which 
translates to a small change (decrease) in a (or equiva- 
lently a s ). However, for higher load levels the value of 
the scale factor a s does change significantly. To see this, 
consider a much higher load level, say F n — 3.5. The 
scale factor for this case is a s is 2.34 x 10 4 , which falls 
well outside the mean value of a s for the limited load 
range F n = 1.446 — 2.26. A plot of the force-displacement 
curve is shown in Fig. [21(a) for F n = 3.5 along with the 
cumulative charge transferred to the substrate £<j. Thus, 
the lack of dependence of a on load is clearly due to the 
limited range of loads studied in Ref. [17]. The plot also 
shows that the model predicts fewer stick-events on an 
average with increase in F n for the same scan length, a 
feature taken as a support for electronic origin of friction 
|17| . The above features emerge due to the separation 
of time scales of the stick (charging) phase and the slip 
(charge transfer) phase. 

As stated earlier, at the charge density level of 
10 8 charges/mm 2 , the slip threshold is controlled by the 



friction threshold A x qtqZ 



3/2 



However, when a is in- 



creased, the onset of the first force drop increases. A plot 
of the force-displacement curve is shown in Fig. [^a) for 
a = 5 x 10 9 charges/mm 2 keeping F n = 2.0. (The cu- 
mulative charge transferred to the substrate E^ is also 
shown.) It is clear that the threshold for the first force 
drop as also the average asymptotic value of F s has in- 
creased from those for a — 10 8 charges/mm 2 . (Compare 
with Fig. ED (a).) 

As stated earlier, the results of the Model II are very 
similar to those of Model I. To see this, we first note that 
the equation for X (Eq. [TOf for Model II differs from that 
of Model I [Eq. ©] by a factor Z 3 ^ 2 . However, the near 
saturation value of Z is almost entirely controlled by Eq. 
©. Thus, the value of r]',, where the instability occurs 
can be estimated by inserting the near saturation value 
of Z in ?7|| = i]'nZ 3 / 2 . This gives a value rf,, ~ 4x 10 6 . The 
rest of the parameters are the same. A plot of the force- 
displacement curve shown in Fig. [3] for F n = 2.0. For 
the sake of comparison, we have also shown the force- 
displacement curve for Model I on the same plot as a 
dashed line. As can be seen, except for a phase factor, 
the two curves follow the same behavior. 

The above calculations are for n = q = 1. However, in 
general, they can be different from unity. Figure 2] shows 
the force-displacement curves for both Model I and II for 
n = 1.75 and q = 2. Again the difference between the 




time 



FIG. 2: (color online) (a) Plots of scaled force and cumulative 
charge transferred to the substrate for F n — 3.5. (b) Plots 
of scaled force and cumulative charge transferred to the sub- 
strate for F n = 2.0 and S m = 0.835 (a m = 5 x 10-4C/m 2 ). 



two models is only in the phase. 
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FIG. 3: (color online) (a) Plots of the scaled force - 
displacement curves for Model II (continuous blue) and Model 
I (dashed red). 

We now consider estimating the value of the scale fac- 
tor a. Since the force and charge transferred curves for 
a given scan length collapse onto a single curve only in 
the asymptotic stick-slip regime [Fig. QJb)], we can use 
x ~ V a valid on an average. Using this in 



/ K\\xdt = a &ddt = a 



{xa t /D)dt (13) 
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FIG. 4: (color online) (a) Plots of the scaled force - 
displacement curves for Model II (continuous blue) and Model 
I (dashed red) for n = 1.75 and q = 2. 



gives 



K\\D ~ airRz asv <j a 



(14) 



where a aS y is the mean value of a (p — > <7 m ) in the asymp- 
totic stick-slip regime. Using z asy and a m leads to a ~ 50 

o 

eV/A for F n = 94 miV. This value is high even after dis- 
counting the ideal nature of the model. 

This prompted us to look for possible contributions 
that affect the scale factor a. A cursory look at con- 
tact charge image in Fig. 2 of Ref. [13] shows that 
the mean charging radius R c is surprisingly large, typ- 
ically 10-20 times the contact radius a with charge den- 
sity <7 ~ 1.6 — 2.4 x 10 8 charges/mm 2 . While the au- 
thors make no comment on such a high value of charging 
radius, one possibility is that charge deposited spreads 
(possibly after the slip event) due to electro-static repul- 
sion [23, 1231 . However, the model uses the natural choice 
that charging area is equal to the contact area for that 
load. Thus, independent of the physical causes for such 
a large charging area, we should either use a larger con- 
tact area j3A n with /3 ~ (-^f) 2 ~ 100 - 400 or use higher 
charge density (by a factor /3) keeping the contact area 
fixed. Using j3A n in place of A n , we get a ~ 0.5 — 0.125 

o 

eV/A. This value is consistent with the reported value. 
Finally, we note that the value of a is nearly the same 
for smooth sliding conditions as x = V a is exact. 

Another factor that affects contact charging is the sur- 
face roughness. The influence of surface roughness on 
friction is well known in tribology. Several models have 
been proposed to calculate the contribution arising from 
roughness [35M37I ] . The topic has been revisited recently 
with a view to accommodate the influence of a hierar- 
chy of roughness scales [38j . However, for the load levels 
considered here, it is easy to show that this is likely to 
contribute about 10% (38[. Here, it may be relevant to 
point out that roughness in this case is likely to be dy- 
namical due to the sliding over long intervals of time [j] . 



V. SUMMARY AND DISCUSSION 

In summary, we have developed a model for the stick- 
slip dynamics observed when the tip of the cantilever of a 
surface force apparatus is dragged on a substrate. Equa- 
tions of motion for the center of the contact area and pen- 
etration depth are set up by considering the tip as a sin- 
gle smooth asperity contacting with a smooth surface for 
which contact mechanics is applicable. Further, the time 
dependent contributions arising from visco-elastic nature 
of the PMMA substrate and possible plastic deformation 
of the softer PMMA are included. The model exhibits 
stick-slip dynamics for a range of values of the param- 
eters. To account for the contact charging at the area 
of contact, these equations are coupled to an equation 
of motion for the total contact charge developed on con- 
tact. Charging occurs during the stick phase and transfer 
of charge to the substrate occurs during slip phase. 

The model captures several features of the experiment. 
First, the correlation between the stick-slip events and 
charge deposited on the PMMA substrate is reproduced. 
This is to be expected as the contribution from contact 
charging ( ~ nN) is orders of magnitude smaller than the 
static friction threshold (~ mN) or even sliding frictional 
threshold. Second, for higher load, the model predicts 
fewer events for the same scan length. Third, lack of de- 
pendence of the scale factor a on normal load observed in 
experiments is traced to the small range of values stud- 
ied in Ref. [17] . Indeed, Eq. (|14p shows that the scale 
factor a is inversely proportional to the area of contact 
for the load. Fourth, in the model, contact charging is 
equal to the contact area for the load. If this relation is 
used, the magnitude of the scale factor turns out to be 
nearly 50 times the experimental value. However, as dis- 
cussed, the experimental charging area is (3 ~ 100 — 400 
times the theoretically calculated contact area for that 
load (see Fig. 2 of Ref. [17J)- Using this value in place of 

the contact area in Eq. (fT4|) gives a ~ 0.5 — 0.125 eV/A. 
This value is close to the measured value. Since the lack 
of dependence of a on load together with the value of a 
has been interpreted to mean that friction is controlled 
by contact charging, it is necessary to check the lack of 
dependence of a over larger range of normal loads before 
the suggestion can be taken seriously. Further, the value 
of a should be nearly same for smooth sliding conditions 
as well. All these features follow naturally due to the sep- 
aration of time scales of stick and slip events. Stick-slip 
events are not affected by contact charging for the charge 
levels in Ref. [l7(. Instead plastic deformation of the 
interface material is the cause of slip. At a mathemati- 
cal level the instability is due to a competition between 
the visco-elastic and plastic deformation time scales with 
the applied time scale. As most experimental features 
are captured, the model provides an alternate explana- 
tion for the results. Lastly, the fact that sliding friction 
is lower than static friction comes out naturally from the 
equations of motion. 

Finally, the model should be applicable to other 
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stick-slip phase observed in frictional sliding experiments 
where plastic deformation is significant, including stick- 
slip observed during scratching of polymer sheets [39M41J J . 



ACKNOWLEDGMENTS 

G. A would like to acknowledge financial support from 
BRNS grant 2007/36/62-BRNS and Raja Ramanna Fel- 
lowship scheme. 



[1] J. Gao, W. D. Luedtke, D. Gourdon, M. Ruths, J. N. 

Israelachvili and Uzi Landman, J. Phys. Chem. B 108, 

3410 (2004). 
[2] A. M. Homola, J. N. Israelachvili, P. M. Mcguiggan and 

M. L. Gee, Wear 136, 65 (1990). 
[3] R. W. Carpick, N. Agrait, D. F. Ogletree, and M. 

Salmeron, Langmuir 12, 3334 (1996). 
[4] A. Berman, S. Steinberg, S. Campbell, A. Ulman and 

Jacob Israelachvili, Tribolgy Letters 4, 43 (1998). 
[5] Y. Mo, K. T. Turner and I. Szlufarska, Nature (London) 

457, 1116 (2009). 
[6] B. Luan and M. O. Robbins, Phys. Rev. E 74, 026111 

(2006); S. Cheng, B. Luan and M. O. Robbins, Phys. 

Rev. E 81, 016102 (2010) and the references there-in. 
[7] K. L. Johnson, K. Kendall and A. D. Roberts, Proc. R. 

Soc. London A 324, 301 (1971). 
[8] K. L. Johnson, Contact Mechanics (Cambridge Univer- 
sity Press, Cambridge, England 1985.) 
[9] I. Szlufarska, M. Chandross, and R. W. Carpick, J. Phys. 

D, Appl. Phys. 41, 123001 (2008). 
[10] F. Heslot, T. Baumberger, B. Perrin, B. Caroli and C. 

Caroli, Phys. Rev. E 49, 4973 (1994); T. Baumberger 

and C. Caroli, Adv. Phys. 55, 279 (2006). 
[11] T. Baumberger, P. Berthoud and C. Caroli, Phys. Rev. 

B 60, 3928 (1999): L. Bureau,T. Baumberger and C. 

Caroli, Euro. Phys. J. E 19, 163 (2006). 
[12] Y. Estrin and Y. Brechet, Pure and Appl. Geophys. 147, 

745 (1996). 
[13] U. Landman, W. D. Luedtke, N. A. Burham and R. J. 

Colton, Science, 248, 454 (1990); U. Landman, W. D. 

Luedtke and E. M. Ringer, Wear 153, 3(1992). 
[14] U. Landman, W. D. Luedtke and J. Gao, Langmuir 12, 

4514 (1996). 
[15] M. R. Sorensen, K. W. Jacobsen and P. Stoltze, Phys. 

Rev. B 53, 2101 (1996). 
[16] C-K Chen, B-H Chen and Y-T Yang, Nanotechnology 

15, 1771 (2004). 
[17] R. Budakian and S. J. Putterman, Phys. Rev. Lett. 85, 

1000 (2000). 
[18] A. D. Roberts, J. Phys D: Appl. Phys. 10, 1801 (1977). 
[19] H. Graf Von Harrach and B. N. Chapman, Thin Solid 

Films 13, 157 (1972). 
[20] A. R. Akande, J. Electrostatics 63, 1025 (2005). 
[21] Lieng-Huang Lee, Journal of Electrostatics 32, 1 (1994). 
[22] J. Lowell and A. C. Rose-Innes, Adv. Phys. 29, 947 

(1980). 
[23] G Ananthakrishna and Jagadish Kumar, in the Proceed- 



ings of II International Workshop on Correlation Microsc- 
tructure - Properties and Multiscale Modelling of Plas- 
ticity (Fuenteheridos, Spain, June 17-21,2009), (Eds.) 
A. Dominguez-Rodriguez, D. Gomez-Garcia, A Munos 
Bernab University of Sevilla, Spain, 2009, pp. 91-94; See 
also Jagadish Kumar, Stick-slip dynamics in peeling of an 
adhesive tape and surface force apparatus, Doctoral The- 
sis, Indian Institute of Science, Bangalore, India, (2010). 

[24] E. Charrault, C. Gauthier, P. Marie, and R. Schirrer, 
Langmuir 25, 5847 (2009). 

[25] D. Maugis, J. Colloid Interface Sci. 150, 243 (1992). 

[26] Z. H. Stachurski, Prog. Polym. Sci. 22, 407 (1997). 

[27] G.Ananthakrishna, Phys. Rep. 440, P 113-259 (2007). 

[28] J. A. Hurtado K-S Kim, Proc. Roy. Soc. London, A 455, 
3363 (1999). 

[29] W. P. Vellinga and C. P. Hendriks, Phys. Rev. E 63, 
066121 (2001). 

[30] A. S. Nowick and B. S. Berry, Anelastic Relaxtion of 
Crystalline Solids New York 1972. 

[31] R. W. Hertzberg, Deformation and Fracture Mechan- 
ics of Engineering Materials, (John Wiley & Sons,Inc, 
New York, 1996). 

[32] G Ananthakrishna and M. S. Bharathi, Phys. Rev. E 70, 
26111 (2004). 

[33] S. Kok, M. S. Bharathi, A. J. Beaudoin, G Ananthakr- 
ishna, C. Fressengeas, L. P. Kubin, and M. Lebyodkin, 
Acta Mater. 51, 3651 (2003). 

[34] A better way of dealing with plastic deformation in the 
x (or in the z) direction would be to write an equation 
of plastic strain rate and couple it to the equation of 
motion for x (or for z). This is not attempted here since 
the equations of motion are already complicated. 

[35] J. A. Greewood and J. B. P Williamson, Proc. Roy. Soc. 
London A 295, 300 (1966) . 

[36] P. R. Nayak, J. Lubr. Technol. 93, 398 (1971). 

[37] A. W. Bush, R. D. Gibson and T. R. Thomas, Wear, 35, 
87 (1975). 

[38] C. Yang and B. N. J. Persson, Phys. Rev. Lett. 100, 
024303 (2008) and the refences therein. 

[39] K. Li, B. Y. Ni and J. C. M. Li, J. Mater. Res. 11, 1574 
(1996). 

[40] W. P. Vellinga and C. P. Hendriks, Tribology Lett. 9, 119 
(2000). 

[41] C. P. Hendriks and W. P. Vellinga, Rev. Sci. Instrum. 
71, 2391 (2000). 



